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of  2  so  that  the  curves  represent  .5,  1.5,  4.5,  and  7.5  dB  of  loss. 

luJ - 

P.N.  MIKHALEVSKY 

\ 


r*! 


TM  No. 
801040 


ABSTRACT 


Detection  threshold  calculations  are  performed  for  cross  correlation 
processing  of  the  received  signals  from  two  fixed  omnidirectional  sensors  of  a 
moving  source  generating  an  arbitrary  random  power  spectrum.  The  effects  of 
doppler  decorrelation  and  volumetric  absorption ,  and  the  degradation  as  a  re¬ 
sult  of  multiple  bottom  bounce  and  surface  reflections  are  included.  The 
equations  for  ROC  curves  are  presented.  A  computer  program  which  accepts  as 
input  the  source  power  spectrum,  noise  power  spectrum,  averaging  time,  fre¬ 
quency  analysis  band,  and  sensor  baseline  separation  solves  for  the  detection 
threshold  given  a  .5  probability  of  detection  and  .001  probability  of  false 
alarm. 
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NOTATION 

Analysis  bandwidth  =  ^2-fl 
Sensor  base  line 
Speed  of  sound  in  water 
Upper  frequency  limit  for  analysis 
Lower  frequency  limit  for  analysis 
Impulse  response  function  accounting  for  volu¬ 
metric  absorption,  surface, and/or  bottom  re¬ 
flections  (s) 

Equals  0  (no  target) ,  or  1  (target  present) 

Noise  at  receiver  i  (i»l,2) 

Transfer  function  of  h(t) 

Noise  spectrum 

Source  function 

Source  power  spectral  density 

Range  from  target  to  midpoint  of  sensor  baseline 
Averaging  time 
Target  3peed 

Variables  of  integration 

Target  speed  along  the  line  of  sound  between 
target  and  receiver  i  (1*1,2) 

Received  signal  at  receiver  1(1*1, 2) 

Volumetric  absorption  loss  coefficient,  function 
of  frequency  and  range 

Time  compression  or  expansion  at  receiver  i 
(1*1,2)  =  1+v^/c 


-5- 


Probability  of  detection 


TM  No. 
801040 


*0 

?„ 

Erf (x) 

c2 

0 


m 


y 


DT 

ML 

SL 

3E 

>T*T 


Probability  of  false  alarm 


A  * 


Error  fraction  3  (2/,Tt) 

Variance  of  the  correlator  output  given 
i  Var{y|k*o} 


variance  of  the  correlator  output  given 
target  is  present  =  Var{y!k*l; 

Expected  value  of  the  correlator  output 

the  target  is  present  3  E(y ; k*l; 

Correlator  output 

Detection  threshold 

Noise  level 

Source  level 

Signal  Excess 

Transmission  Loss 


no  target 


the 


given 
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INTRODUCTION 

Detection  threshold  calculations  are  performed  for  cross  correlation 
processing  of  two  received  signals  each  having  propagated  via  no,  one  or 
many  bottom  bounces  and/or  surface  reflections  from  a  moving  source  genera¬ 
ting  an  arbitrary  random  power  spectrum.  The  two  receiving  sensors  are  omnidirectional, 
fixed  and  of  arbitrary  baseline  separation,  a  binary  hypothesis  test  is 
performed  and  the  detection  threshold  computed  for  a  probability  of  de¬ 
tection  of  .5  and  a  probability  of  false  alarm  of  .001. 

The  effects  of  doppler  and  absorption  on  the  mean  correlator  output 
are  investigated.  It  is  shown  that  in  general  degradation  due  to  doppler 
will  be  >2d3  when  3Tj  (c/V)  (R/b)  { 1+1/4  A)2}*3.  This  can  only  be  considered 

a  rough  guide,  however,  as  the  doppler  degradation  is  a  function  not  only  of 
analysis  bandwidth  but  upper  and  lower  cut  off  frequencies.  The  detection 
threshold  curves  in  fact  exhibit  high  peaks  corresponding  to  nulls  in  the 
correlator  output  as  a  function  of  range  as  a  result  of  this  dependence. 

For  an  analysis  band  of  300  Hz  it  is  shown  that  absorption  is  <3dB  for  ranges 
<30km  and/or  frequencies  <600  Hz. 

Bottom  losses  axe  taken  from  measured  values  in  the  abyssal  plains 
south  of  3ermuda  and  extrapolation  based  upon  other  measurements.  Surface 
interactions  are  characterized  by  the  Kirchoff  approximation  assuming  a 
random  Gaussian  surface,  though  computations  of  detection  threshold  for 
surface  bounce  paths  are  not  performed.  The  purpose  of  this  report,  however, 
is  not  to  characterize  surface  and  bottom  interactions  but  how  to  account  for 
their  effects  in  cross  correlation  processing  given  specific  conditions  and 
assumptions . 
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ANALYSIS 

In  general  we  have  two  received  signals 

=  ks£*rfc)*  -»-«,&)  a) 

(fc)  *  k s  C«iO  +  n-J-O  *wt±)  (2) 

where  s(t)  is  the  source  function,  n. (t)  is  the  noise  at  receiver  i,  h(t)  is 

i 

the  impulse  response  function  accounting  for  volumetric  absorption,  surface  and/ 
or  bottom  reflection(s) ,  w(t)  is  a  weighting  function  corresponding  to  a  band 
pass  filter  in  the  frequency  domain,  and  k  »  0  when  no  target  is  present  and 
k  ■  1  where  the  target  is  present.  The  *  denotes  convolution. 

<*  ,  =  1+  V,  /C  Amd  ,  ' 3> 

<*a  =  1  (4) 

where  v  is  the  target  speed  along  the  line  of  sound  at  sensor  i.  The  effect 
of  doppler  is  thus  a  time  compression  or  expansion  as  indicated  in  Eqs.  (1) 
and  (2).  It  is  also  assumed  that  s(t),  n^(t),  and  n2<t)  are  zero  mean,  un¬ 
correlated,  Gaussian,  and  stationary. 

The  maximum  output  of  the  correlator  occurs  when  x_^  is  shifted  in  time  by 
the  appropriate  time  delay  relative  to  x^  which  assuming  this  has  been  ac¬ 
complished  we  can  write  as,  ^ 

y  9  y  J  *■£  (5) 

'lx 

For  T3  large,  independent  of  the  distributions  of  x^  and  y  will  be 
a  Gaussian  random  variable,  thus  knowledge  of  the  mean  and  variance  of  y  is 
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all  that  is  required  to  completely  characterize  the  probabilistic  behavior 
of  y. 

A.  The  Expected  Value  of  y 

The  expected  value  of  y  given  k  recalling  the  above  assumotions  is 

E  [y i  K]  3  J-  $  k1  E  [[  S^Vt)  *  Uft )][  bUjz)  *vr(t)*  d±  { b ' 

If  we  perform  the  indicated  convolutions  and  write  the  correlation  function 
cf  s(t)  as  the  inverse  of  its  power  spectral  density  G  (f)  we  obtain, 

T/z  *«  =°  oo 3 


c[yik]  =  7  j  *.t  j  ay,  { ci*2  [  «y3  $  <xy4  J  a-f 

— _  co  -co  -co  _  oo  „  oo 


OO  —  OO  _  oo  .00 

l  2TT -j  (oC  -  *^1^1  “  y  I  ^3  ”  +^2.^2 ) 


(7) 


•  C-,«)  a 

•  w'x.'w-fxO'h^JkCXfl) 


Integrating  over  the  x.'  s  we  have 

r;/*  r°°  izrr-f  (*,t  - 

j  )  ci-f  GSC4)  e 

T/z  -<» 


(8) 


where  H(f)  is  the  transfer  function  accounting  for  bottom  and  surface  inter¬ 
actions  and  volumetric  absorption  and  w(f)  is  a  band  pass  filter  with  a  gain 
of  1/  J~2B  for  -f2<f<-f1  and  f^<f<f2  and  zero  elsewhere.  Interchanging  the 

order  of  integration  and  integrating  over  t  we  have , 

co 

&:yik]=  [tt£t  («W2)]  Gs  '-f) 

-  00 

'  w  ^r?)w  4) 


(9) 
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we  can  write , 
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is. 

$  S,ff)  SMC  [rr-frU 

f. 


|HCP|  df 


CO) 


,  s±n(x) 

where  sinc(x)*  — - -  . 

This  is  as  far  as  we  can  proceed  without  knowledge  of  G^Cf)  and  H(f) . 

One  can  readily  see  from  Sq.  (10)  that  the  effect  of  doppler  is  a  sine  function 
weighting  of  the  source  spectrum.  The  sine  function  has  its  first  zero  crossing 
at 


i  ■  Lt(1i -.<*■>] 


an 


If  the  center  of  the  analysis  band  is  near  or  greater  than  fQ  and  B>fQ  severe 

degradation  (>>3DB)  of  the  correlator  output  can  be  expected.  If  we  choose 

a  "worst  case"  geometry,  namely  the  target  is  travelling  an  a  course  parallel 

to  the  sensor  baseline  and  crossing  the  perpendicular  bisector  of  the  baseline, 

then  the  condition  3>f  can  be  written  as 

c 

BT  +  <UI 

As  a  representative  case  the  loss  in  the  mean  due  to  doppler  is  plotted 
in  Fig.  1,  assuming  the  source  spectrum  is  flat  and  the  analysis  band  extends 
from  0  to  B  Hz.  Eq.  (10)  is  used  neglecting  the  bottom  transfer  function. 

For  bottom  or  surface  bounce  paths  doppler  degradation  will  not  be  as 
severe  as  for  direct  path  propagation  as  the  target  speed  V  is  reduced  by  the 
factor  cos6  where  5  is  the  bottom  or  surface  grazing  angle.  Furthermore  though 
it  is  true  that  the  doppler  degradation  will  be  more  severe  as  the  range  is 
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reduced,  the  signal  is  experiencing  less  transmission  loss  which  cam  reduce 
or  negate  this  effect.  More  discussion  on  this  will  be  given  in  Section  E. 
However,  because  the  sine  function  does  go  negative  it  is  possible  to  obtain 
zero  output  if  the  analysis  band  is  appropriately  placed  relative  to  the  main 
and/or  side  lobes  of  the  sine  function.  Because  the  zero  crossings  of  the 
sine  function  will  change  as  the  range  to  the  target  changes  (relative  doppler 
between  sensors  is  a  function  of  range)  nulls  in  mean  correlator  output  and 
therefore  corresponding  spikes  in  detection  threshold  appear  as  functions  of 
range  ( see  Section  S . ) 

Assuming  no  surface  and  bottom  interactions  we  can  express  Hif)  as 


r~  ^  _  *1 

y^Ci)  =  Jixp  j 


(13) 


where  a^(f,R)  is  the  absorption  loss  coefficient  which  is  a  function  of  frequency 
and  range  (see  for  example  ref.  (!}).  If  we  apply  Sc.  (13)  to  Sa.  (1C), 


assume  a  flat  source  spectrum  and  no  doppler  we  have, 


(14) 


*1 


for  the  Mean  Absorption  Loss.  Note  that  this  quantity  is  range  dependent. 
Eq.  (14)  is  plotted  in  Fig.  2  with  B®300  Hz  for  low  frequency  cutoff  f^  vs. 
range  B. 

3 .  The  Variance  of  y 

2 

The  expected  value  of  y  given  k  is 


'Z  .2. 

■  [y"  1*^1  s  Za.  S  5  ^  [*R.i  )  ^Ri ' -i)  2.)  J  G.'tj 


(IS) 


’T/l  ^ 
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Exploiting  the  Gaussian  assumption  we  moment  factor  Eq.  (15) , 
'4.  '/2. 

E f  y* !  '< j  *  ZL  j  ^  [  -  [**.» ; ^ L^ri  ^2} 

T  -t/2-t4 


(16) 


+  £  L^ftj  &0  ^Rl  ^2^1  E  E*Ra^  * 


Thus  the  variance  of  7  is 

T4  ,TA 


(17) 


Vav-Tyik]  =  -  ^  j" 

-r/z  -T4  (1 

*+•  c  [%R|  Cz ,)  4Xfti^r*^]  £  E^R2  )jj  ctt  1  cfcfcj. 

Applying  Eqs.  (1)  and  (2)  in  Eq.  (17),  performing  the  indicated  convolutions, 
transforms,  and  changing  the  double  integration  over  t^  and  t  to  a  single 
integration  over  Tst^~t2  w^ere  possible  we  have, 


Vw  V  k  j  -  i  f  0  -  7 )  {  ®®  -©®+€®-®©} 

' T  T4  T4 


(18) 


$  ©@  cfe;ctz 


where 


,x  A 


~Vz  ~T/- 


®  '  "b"  ^  G-J^)  Ccri,  (2TT$ oi{r) 

®  =  T  i  G*s  (4)  Coa,  (zrr^^r)  1KCP\”  M 

,  ^ 

0  -  ^  j  e5(f)  CBt[2nf 
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A 

©  -  T  j  Gs  ff)  cos  r^TT  4  df 

f, 

®  -  i)  Sv,  C4)  cos  ferfr)  .If 

4,  ' 

®  *  S  $  G*  &)  C05  ^2Tt4t)  <4 

ft  7- 

The  approximations  arise  as  a  result  of  the  assumptions  that  |  "^x~c  t  ^  ® 

and  lf.fi  <‘3. 

C.  Binary  Hypotheses  Test 

We  recall  that  for  TB»1.  the  correlator  output  y  is  a  Gaussian  random 
process  when  the  target  is  present  (k»l)  or  not  (k=0) .  Applying  a  binary 
hypothesis  test  we  obtain 


_  1  “I  .  -  (^1  1 

7-n  0-  zPo)  =  ^  Er* 


(19) 


where  Erf  ^(x)  is  the  inverse  error  function,  j2  »  Var{y|k*0)  ,o2  ■  Var{y!k*lK 

0  I 

5  S{y | k*l: ,  PQ  is  the  probability  of  detection, and  P?  is  the  probability 
of  false  alarm.  Eqs.  (10),  (18),  and  (19)  define  the  Receiver  Operating 
Characteristics  for  cross  correlation  processing,  with  PD*.S,  and  P?*.001 
we  obtain  from  Eq.  (19) , 

i,»  ra-lz  £-rf(.w)  =  <r.v'3.od|...) 


(20) 


m, 


where. 


0*5  C4) si'ncjrr^Tfoc,-^ )][  J >T(-f)|  cL-f 

4, 


(21) 


where  G  (f)  is  the  source  expectrum  attenuated  equally  by  frequency  such  that 
s 
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Eq.  (20)  is  satisfied,  and 

A iA 


-  g*T  I  ^  ~ T*  ^  i"  Gnjfd  CAd^rrfr)  GsS  (2TT  vr)  aif  elvolr 


(22) 


If  we  assume  equal  noise  at  the  receivers  and  exploit  the  fact  that  TB»1, 


<Ti 


i 


*  ao 

WT  J  I  W«H 


W 


G> 


(23) 


where  HMf)  is  a  passband  filter  of  gain  1  for  '  fj_  '<£< '  f2i  •  Noting  that  the 
inner  integral  is  the  transform  of  R  (t),  the  correlation  function  of  the 


noise,  convolved  with  h  (r) ,  and  then  applying  Parseval's  theorem  we  have 


simply, 


[&*&)]  dj- 


(24) 


If  we  assume  a  flat  signal  (Gg(f)=S)  and  noise  (G  (f)»N)  spectrums ,  no 
doppler ,  bottom,  surface,  or  volumetric  losses  and  combine  Eqs. (20) ,  (21)  and 
( 24 )  we  have 


^SN  »  Y==  C.^8) 


(25) 


where  R^is  the  signal  to  noise  ratio  (;S/N) .  This  result  was  first  obtained 
by  Nuttall  r2}. 

0 .  The  Transfer  Function  3(f) 

In  general  for  multiple  bottom  and  surface  reflections  we  have 


for  N  bottom  and  M  surface  reflections.  The  exponential  term  accounts  for 
volumetric  absorption  and  R^,  is  the  total  range  including  the  bounces.  Fur¬ 
thermore,  the  Hg's  and  H^'s  are  in  reality  average  filter  responses,  averaged 
over  the  time  and/or  specially  varying  surface  as  the  bottom  and  sea  surfaces 
are  random  surfaces.  These  realities  in  fact  limit  the  applicability  of  the 
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of  the  analysis  in  Section  3,  particularly  the  moment  factoring  used  to  derive 
the  variance,  which  requires  essentially  that  the  result  of  the  convolution  of 
s(t),  a  Gaussian  random  process,  with  h(t)  the  impulse  response,  also  be  a 
Gaussian  random  process.  Note  that  independent  of  the  nature  of  H(f) , 
the  correlator  output  will  always  be  Gaussian  for  TB>>1. 

If  we  assume  the  reflecting  surface  is  random  with  displacement  I(x,t) 


and  employ  the  ICirchoff  approximation  we  cm  write  (3}, 

x3.Cf)»  5  p  cy  ex?  e* )  <*  tj 


(27) 


where  p_  is  the  probability  density  function  of  the  surface  for  the  jth 

J 

reflection,  $  is  the  speed  of  sound,  and  i ^  is  the  surface  grazing  angle. 
Assuming  a  Gaussian  surface  we  have, 


K*Cf>*  «P 


(29) 


where  a,2  is  the  variance  of  the  surface  for  the  jth  reflection.  With  this 
characterization,  extension  to  multiple  reflections  is  straightforward. 

For  a  single  surface  reflection 

£  [yIk]  s  b  J  Gift)  si»e  (rfT  (a,-  «z)]  (29) 

“f" » 

•  **p  {-  <xA  (Wtf  ,R )  -  *  frrf.11)  -  2  [  fr?+ai  )]}  M 

where  it  i3  assumed  that  a.2  and  9  are  the  same  for  the  surface  reflection  to 

;  s 

each  of  the  two  receivers.  If  we  make  the  further  assumption  that  the  two  re¬ 
flection  surfaces  are  separated  sufficiently  such  that  they  are  statistically 
uncorrelated  (with  the  same  variances)  computation  of  the  variance  of  the 
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correlator  outputs  is  straightforward  with  similar  expressions  for  H(o.f)H*(a2f) 
in  Eq.  (18)  as  given  in  Sq.  (29) . 

A  similar  approach  can  also  be  applied  for  bottom  reflections.  However, 
when  this  method  was  applied  predicted  and  observed  bottom  losses  did  not 
agree,  with  less  bottom  loss  predicted  at  low  frequencies  and  more  bottom  loss 
predicted  at  ^igh  frequencies  than  observed.  Using  measured  values  {4}  for  an 
area  south  east  of  Bermuda  in  which  a  prototype  of  the  cross  correlation  pro¬ 
cessor  is  operating  single  bounce  losses  were  extrapolated  based  on  these 
and  other  bottom  loss  measurements  (1).  The  curves  are  shown  in  Fig.  3  for 
H3(f)!2  vs.  grazing  angle  at  various  frequencies.  Clearly  surface  and  bottom 
interactions  must  be  carefully  investigated  for  the  area  of  operation  or  else 
grievous  errors  could  be  made  in  predicting  system  performance. 

2.  Calculation  of  Detection  Threshold 

A  computer  program  has  been  developed  which  accepts  as  input  the  target 
spectrum,  target  speed,  upper  and  lower  frequency  limits  for  the  analysis 
band,  and  the  averaging  time,  and  computes  the  detection  threshold  vs.  range. 

The  noise  spectra  are  taken  from  Wenz  {5}  with  three  spectra  available: 
heavy  noise  (shipping  density  1,  sea  state  5) ,  average  noise  (shipping  density 
.45,  sea  state  3),  and  light  noise  (shipping  density  0,  sea  state  0).  Eqs. 

(20) ,  (21) ,  and  (24)  are  used  with  numerical  integration  employing  the 
Romberg  Quadrature  technique  {6}  to  compute  the  detection  threshold, 

DT  5  SL  -NL  (30) 

o 

where  SL  is  the  signal  level  which  satisfies  Eq.  (20) ,  and  NL  is  the  noise  level, 
o 

For  broadband  signals  I  define  these  terms  as  follows: 

SL  »  10  log, (31) 

O  10  1 
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GRAZING  ANGLE  DEGREES 

Fig.  3.  Bottom  loss,  |tl  (f)|2,  vs.  grazing  angle. 
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and 


NL 


llGn(f)df 


(32) 


Similarly  the  target 


source  level  is  defined 

G  (f)df 


(33) 


In  Figs.  4-7  are  plotted  examples  of  detection  threshold  and  signal  excess 
calculations  for  the  double  HX-29  {7}  source  towed  at  IS  kts.  in  average  noise 
conditions  assuming  a  sensor  separation  of  1.5  km.  Signal  excess  is  computed 
assuming  a  spherical  spreading  loss , 


3E  ■  SL-NT.-DT-20  log^R  (34) 

The  source  spectrum  is  given  in  Fig.  8.  This  is  the  maximum  source  level  for 
a  driving  bandwidth  of  one  Hz  swept  across  the  frequency  band.  For  broadband 
operation  the  source  spectrum  shape  remains  the  same  but  must  be  reduced  by 
10  log10B‘  The  analirsis  band  ^  Figs.  4-7  is  150  Hz  to  300  Hz  with  an 

averaging  time  of  2.56  secs.  Figs.  4  and  5  are  the  detection  threshold  and 
signal  excess,  respectively,  for  direct  path  propagation.  Figs.  6  and  7  are 
the  detection  threshold  and  signal  excess  respectively,  for  a  one  bottom 
bounce  path. 

Note  the  detection  threshold  peak  in  Fig.  4  which  is  the  result  of  a 
corresponding  correlator  output  null  due  to  doppler.  For  larger  analysis 
bands  and  longer  averaging  times  there  can  be  many  such  detection  threshold 


assuming  direct  path  propagation. 
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Pig.  5.  Signal  excess  for  the  case  shown  in  Pig. 
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RANGE 

excess  for  the  case  shown  in  Fig. 
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of  the  double  HX-2?  gource  driven  in  &lHz  band  across 
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peaks.  However,  for  the  bottom  bounce  path  (Fig.  6)  the  cosine  reduced 
velocity  along  the  line  of  sound  eliminates  this  peak  (see  Section  A) . 

Comparison  of  Fig.  4  with  Fig.  5  and  Fig.  6  with  Fig.  7  reveals  that 
the  reduction  of  transmission  loss  as  a  result  of  reducing  the  range  more 
than  compensates  for  the  loss  due  to  doppler  at  the  shorter  ranges  (except 
of  course  for  the  null.) 


